Thermal transport in a Luttinger liquid 
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We study thermal transport in a one-dimensional (Id) interacting electron gas, employing the 
Luttinger liquid model. Both thermal conductance and thermopower are analyzed for a pure Id gas 
and with impurities. The universal ratio of electrical to thermal conductance in a Fermi-liquid - 
the Wiedeman-Franz law - is modified, whereas the thermopower is still linear in temperature. For 
a single impurity the Lorenz number is given by L{T 0) = 3Lo/{2g + g^) - with Lq the Fermi 
liquid value - and the conductance 1/2 < g < 1. For g < 1/2 the Lorenz number diverges as T — > 0. 
Possible relevance to thermal transport in conducting polymer systems is discussed. 



PACS numbers; 72.15.Jf 71.27.4-a 

The Wiedemann-Franz law, which relates the thermal 
and electrical conductivity {k,<7) of metals, played a cen- 
tral role in the historical development of the quantum 
theory of solids. The Lorenz number, L — k/uT, origi- 
nally computed within classical Drude theory, gave fortu- 
itous agreement with experiment due to cancelling errors. 
The quantum theory corrected the errors, and improved 
the agreement. For non interacting electrons Chester et. 
al. [|] showed that the Lorenz number is given exactly 
by Lq — (7r^/3)(A:B/e)^, for arbitrary impurity scatter- 
ing strength. In the 1980's Castellani et. al. |^ argued 
that this universal value was robust even with inclusion 
of electron interactions, provided the system remained 
metallic. Thus, a universal Lorenz number appears to be 
a defining characteristic of the Fermi liquid phase. 

In recent years there has been tremendous interest 
in conducting phases which are not Fermi liquids. A 
paradigm for these are Id interacting electron gas mod- 
els, which exhibit a non Fermi liquid phase even for weak 
interactions ||]. The resulting Luttinger liquid phase is 
characterized by a dimensionless conductance, g, which 
controls various power laws, such as the singularity in the 
momentum distribution function. The resurgence of in- 
terest in the Id Luttinger liquids stems both from the re- 
cent ability to lithograghically pattern true one-channel 
quantum wires 0| |^ and from the realization that Id 
edge states in the fractional quantum Hall effect are Lut- 
tinger liquids Other non- Fermi liquid phases arise 
in quantum impurity problems |^], such as the multi- 
channel Kondo model which is possibly relevant to heavy 
Fermion materials. Bulk 2d non-Fermi liquid phases have 
also been suggested in compressible Hall fluid phases ||] 
and in the cuprate superconductors Q. 

It is natural to anticipate that thermal transport in 
such non-Fermi liquid phases will be qualitatively differ- 
ent, and might help characterize and distinguish them 
experimentally. In this paper, we consider in detail ther- 
mal transport in the Id Luttinger liquids. We show that 
the Lorenz number can be substantially modified from 
its Fermi liquid value, Lq. The thermopower, Q, on the 



other hand, shows characteristically metallic behavior, 
Q — cT. As in conventional metals, the coefficient c is 
non-universal, depending on the curvature of the energy 
bands and the energy dependence of the scattering rates. 

While thermal transport measurements in quantum 
wires and quantum Hall samples are undoubtedly ex- 
tremely challenging, a remarkable recent experiment has 
demonstrated the feasibility of such experiments [ p^ . 
Thermal transport measurements in bulk quasi Id sam- 
ples, such as conducting polymers, are much easier, but 
3d crossovers may tend to complicate the analysis. 

Pure Luttinger liquid: We begin with a model for an 
interacting spinless Id electron gas in the absence of any 
impurities, which has a bosonizcd Hamiltonian density 



Ho = nvo{Nl + Nl+ 2\N+N^). 



(1) 



The right and left moving electron densities, N±, satisfy 
Kac-Moody commutation relations: 



[N±{x),N±{x')] = ±{i/2n)dj{x~x'). 



(2) 



The interaction term mixes right and left movers, but 
can be shifted away as usual by defining new fields 



N± - [g{n+ + n_) ± {n+ - n^)]/2g 



(3) 



with A = (1 — (7^)/(l + g^)- In terms of n± the Hamilto- 
nian decouples into right and left moving sectors: 



with renormalized velocity, v = {2gvo) / {1 + g'^) 
fields also satisfy a Kac-Moody algebra, 

[n±{x),n±{x')] = ±{ig/2TT)d.J{x - x'). 



(4) 

The new 
(5) 



Consider now transport in an ideal Luttinger liquid. 
Intitally, we ignore additional anharmonic interaction 
terms (eg. n^rt_) which couple the right and left mov- 
ing modes in (4). (For chiral quantum Hall edge states, 
these will be absent.) It then suffices to consider only 
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a single right moving channel, n = n+. Such an ideal 
chiral channel can be characterized by transport coef- 
ficients which relate changes in the electrical and 
thermal currents to changes in the chemical potential, 
and temperature T. These coefficients are equivalent to 
Landauer two- terminal transport coefficients , defined 
with "ideal" reservoirs. For ideal quantum wires, these 
coefficients are not measured directly, since the contacts 
do not couple selectively to right and left moving modes. 
However, they can be measured directly for quantum Hall 
edge states. 

The charge density in a chiral channel is conserved by 
(4), and satisfies dtn + dxJ = with an electrical current 
J = vn. Changing the chemical potential, fi, alters the 
electrical current. Balancing the energy in (4) with 
a —fin term, gives AJ = {g/2T:)Afi, or upon restoring 
units an electrical conductance, G — C^^ = ge^ /h. 

Heat carried by a chiral channel is likewise conserved 
by (4). The continuity equation dtug + dxJq = is 
satisfied by the thermal energy density uq — {Trv/g)n'^ 
and thermal current Jq = vnq. The thermal energy at 
temperature T can be expressed in terms of the chiral 
Luttinger modes as 

f°° 

n-Q = / TT^kbuj^ (6) 
Jo 

with — (e*^" — 1)^-^ and cok = vk. This gives 
Jq = (7r^/6)(fc_BT)^//i, and leads to a "quantized" ther- 
mal conductance, K = C"^^ = dJg/dT = [n'^ l'i)k%T /h. 

For an ideal Luttinger liquid we can then define a "two- 
terminal" Lorenz number 

Adcai - K/TG = L0/5. (7) 

For g = 1, we recover the Fermi-liquid value Lq = 
(7r^/3)(fc_B/e)^. With repulsive interactions {g < 1) the 
Lorenz number is larger. 

The off-diagonal transport coefficient C^^ = dJ/dT, 
which determines the thcrmopower, is zero within the 
present model, due to the implicit linearization of the 
electronic bandstructure near the Fermi energy. The 
effects of dispersion can be included via the third or- 
der interaction term, 7ii„t — An^, which is normally ig- 
nored because it is formally "irrelevant" . The coefficient 
A is proportional to the change in Fermi velocity with 
chemical potential, dv/dfi. The resulting thermopower, 
Q = C^'^ /C^^, is linear in temperature, 

Q = -{7r^kl/3gev){dv/dti)T. (8) 

In quantum wires, anharmonic interactions ignored 
above will couple the right and left moving modes. The 
right and left moving thermal currents will no longer be 
independently conserved. However, in a translationally 
invariant system, thermal currents cannot fully relax due 
to constraints of momentum conservation. Such anhar- 
monic interactions might nevertheless effect the value of 
K. Umklapp processes would allow a decay of thermal 
current, but freeze out at low temperatures. In any event, 



impurity backscattering will dominate these interaction 
effects in the thermal resistance of real quantum wires. 

Single Impurity: We now consider a single impurity in 
an otherwise ideal Luttinger liquid, as a first step toward 
inclusion of many impurities. (A single impurity is also 
relevant to point contact experiments in the quantum 
Hall effect.) A weak potential scatterer at the origin can 
be modelled by adding a term to the Hamiltonian, 

T-Csack = -ts cos(0+ - cf)^)5{x), (9) 

where ts is the amplitude for 2kp electron backscatter- 
ing. This process has been expressed in terms of the bo- 
son fields (/)±, related to the densities n± — ±dx4>±/'2,T:. 
The operator exp(i0_|_) creates an excitation with frac- 
tional charge, ge. Thus each backscattering process re- 
flects fractional electron charge. 

An impurity which strongly backscatters can alterna- 
tively be modelled as a tunnel junction between two de- 
coupled semi-infinite Luttinger liquids In this case, 
the chiral density n+ can be taken to describe the right 
and left moving pieces of one semi-infinite Luttinger liq- 
uid. The appropriate term which tunnels an electron 
(charge e) through the junction is then, 

Utunn = -t COS. - (j)-) I g) 6 {x). (10) 

To proceed we follow Fendley et. al. jlj] and first de- 
fine new fields which propogate in the same direction: 
<j)i{x) = 4>+{^) ^'^d 4'2{x) = (/)_(— a;), and associated den- 
sities, rij = dx4ij with j — 1, 2. One can then define com- 
muting even and odd densities, n = n+ — n_ = dx'Pl'^'''^ 
and N = n+ -I- n_ . The full Hamiltonian with backscat- 
tering present factorizes, 

n = {-KV /2g){n'^ + N"^) - tB cos(j)5{x). (11) 

The backscattered electrical current is given by 

J= dtn/2 = gtB sin (l){x^O), (12) 

J X 

where the second equality follows from commuting n with 
the Hamiltonian. For the case of a tunnel junction, the 
tunnel current is J = tsin((j){x = 0)/g). Similarly, the 
backscattered thermal current can be written, 

Jq= f dtiK - ^o")/2 = i7Tv/g)N{x = 0)J, (13) 

where again the time derivatives are evaluated by com- 
muting with the Hamiltonian. This form also holds for 
the tunnel junction. Notice that the thermal current has 
been decomposed into a product of two commuting con- 
tributions: The even density TV, and the electrical current 
J which depends only on the odd boson. This remarkable 
simplification, enables us to derive an expression relating 
the thermal and electrical conductances. 

To this end consider the current correlation function, 

p^it) = iQify < [j{t), j(o)] >= ieit)iP>it) - p<{t)) 

(14) 
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from which the electrical conductance follows: 

ReG{u;) = -ImP^iuj) = ^{P>{lj) - P<(t^)). (15) 

Denoting the corresponding correlators for the thermal 
current, Jq, with a subscript Q, the thermal conductance 
can be obtained from: 



K = lim —ImP^(uj) = -^P;^(uj = 0), 



(16) 



where the latter equality follows upon using the detailed 
balance relation, -Pq (w) — exp(/3Li;)Pg (w). 

The relation between the thermal and electrical cur- 
rent operators (13), allows us to relate their respective 
correlation functions: 



P<{t)^{nv/grD<{t)P<{t) 



(17) 



where D2 is the even boson density-density {N — N) cor- 
relation function. Using the fact, D2{t) — D^(— i), al- 
lows one to express the thermal conductance as: 



l.UTT.o f du) , , „ , , , 



(18) 



The function D2 can be readily extracted since Ti in 
(11) is quadratic in N, giving D2 (uj) — {g/nv'^)uje^^buj- 
In addition using (15), along with detailed balance, 
relates P^ to the electrical conductance: P'^{uo) — 
2ojb^ReG{u!). We thereby obtain our final expression 
relating the thermal and electrical conductances through 
the impurity: 



K = 



1 



8.gT2 



dco 



sinh^(w/2T)' 



(19) 



In the absence of any backscattering, ReG{uj) = g/2TT, 
which gives the pure result K — {tt'^ /3)kgT/h. In the 
limit of strong backscattering, K can be obtained from 
(18) by calculating P'^(w) perturbatively in powers of 
t, the electron amplitude to tunnel through the junc- 
tion. To leading order, P is the correlation function of 
the current, J — tsm((j)/g), evaluated with the free odd 
boson Hamiltonian (ie. (11) with = 0). One finds 
c a non-universal constant depend- 
ing on a short-time cut off. This same constant also en- 
ters the electrical conductance, and so drops out in the 
Lorenz ratio, which in this limit is found to be: 



L = 3Lo/{2g + g^). 



(20) 



For non-interacting electrons {g = 1) this reduces to the 
Fermi-liquid value, but with repulsive interactions {g < 
1), is larger. 

As one lowers the temperature, the Lorenz number 
crosses over from the pure value, L = L^/g to the strong 
backscattering value (20). For the special case oig — 1/2, 
an expression for this crossover can be obtained explicitly. 
In particular, when g = 1/2, a closed form expression for 
the a.c. conductance follows from the exact solution [^: 



ReG{uj) 



2ttuj 



dE{fE - fE+u.) 



(21) 



with = {e^^ + 1) ^. Here Tb is a crossover temper- 
ature scale, Tb ~ 4''^^^^^ ~ *b (for = 1/2)- Together 
(19) and (21) allow one to compute the Lorenz number 
for arbitrary T/Tb- For T >> Tb one finds L = 2Lq in 
agreement with the pure result (7), but for T << Tb, 
the result is L = 18-Lo/5 - a factor of 3/2 larger than the 
strong backscattering result (20). This discrepancy can 
be traced to an irrelevant operator ignored in the pertur- 
bative calculation leading to (20), but included implicitly 
in the 5 = 1/2 crossover. Specifically, consider a pertur- 
bation coupling the electron densities across the junction 
in the strong backscattering limit: 



(22) 



This term does not transfer charge across the junction, 
but does transfer energy and so contributes to thermal 
conduction. Moreover, it feeds into the a.c. electrical 



conductance [g^ as ReG{ijj) 



Insertion in (19) 



then gives a contribution to the thermal conductance 
varying as if a^T^, which must be added to the T^/g-i 
term coming from electron tunnelling. 




12/5+C 



3/(2g+f ) 



FIG. 1. Lorenz number L for transport though a single 
impurity versus temperature, for several values of g. 

Three cases should then be distinguished, as sketched 
in the Figure. For g > 1/2 the electron tunnelling term 
dominates the thermal transport at low temperatures, 
and (20) is correct as T — > 0. For the soluble case g = 1/2 
both processes vary as T^ and contribute to the Lorenz 
number as T ^ 0. Since the coefficient "a" depends on 
details of the junction, a non-universal Lorenz number is 
predicted in this case. However, since the contribution to 
K from (22) is positive (proportional to a^), the Lorenz 
number should be bounded below by (20): K = 12/5 + G 
with C > 0, consistent with the exact solution. Finally, 
for g < 1/2, the interaction term (22) dominates K at low 
temperatures and the Lorenz number diverges as T — > 0: 



{a^T^)/T'^3-i ^ 



(23) 



Physically, the inclusion of strong electron interactions 
enables heat to be transported across the junction much 
more readily than charge. 
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To compute the thermopower associated with trans- 
port through the impurity it is necessary to include in 
the Hamiltonian (11) terms which break particle/holc 
symmetry. In addition to bulk cubic interactions aris- 
ing from dispersion, another local term is of the form 
Ncos{(f>)6{x), which arises from an energy dependence of 
the matrix element Ib- One finds a thermopower linear 
in temperature, Q = cT, with non-universal coefficient. 

The present results can readily be generalized to in- 
clude the electron spin degree of freedom. Assuming 
SU(2) spin symmetry, the Luttinger liquid can be charac- 
terized by a dimensionless charge conductance t/p, which 
is equal to 2 for non-interacting electrons (5|. For Qp > 
2/3 the Lorenz number L is found to cross over from 
io(2/5p) to 3Lo(ffp + 2)^/(8gp(gp + 1)) as one lowers the 
temperature and scales from the weak to strong backscat- 
tering regimes. For gp < 2/3, L diverges as T^~^/9p in 
the low temperature limit. 

Many Impurities: When many impurities are present 
in a repulsively interacting Id Luttinger liquid, the con- 
ducting state is unstable to localization. However, in sys- 
tems with sufficiently dilute but strong scatterers, there 
will be a range of temperatures over which the above 
single impurity results should be observable. Specifi- 
cally, consider a model of many Id conductors, which 
are coupled electrically by dilute weak- links, with a typ- 
ical large separation L. These weak-links would serve as 
"impurities" in the Id transport. Such a model might 
be appropriate for some conducting polymer systems, in 
which the polymer backbones provides the Id conduc- 
tion channel, and the weak-links arise from tunnelling 
between different polymers at their ends. Applicability 
of this model requires the temperature to be above the 
tunnelling rate, tj_, between parallel neighboring poly- 
mers. Otherwise the transport would no longer be Id, but 
would take place coherently across many chains. More- 
over, we require that tunnelling events across successive 
weak-links be incoherent, so that ksT > fiv/L, where 
V is the Id Fermi velocity. In this temperature regime, 
the bulk electrical conductivity should vary as a power 
law, ct(T) - r2/9-2, before crossing over at lower tem- 
peratures either into a localized regime or a bulk metallic 
state. Our single impurity results imply that for g > 1/2, 
the thermal conductivity should also vary as a power law, 
T^/f-i, with a universal non-Fermi liquid Lorenz 
number given by (20). For g < 1/2 an even larger (and 
diverging) Lorenz ratio is predicted. 

In some conducting polymer samples, the electrical 
conductivity does vary as a power law with tempera- 
ture iQ. This power law has been interpreted ||l^ as 
being in the vicinity of a bulk metal-insulator transition. 
However, at a 3d Anderson localization transition a tem- 
perature independent thermopower is predicted JTst , in 
contrast to the measured behaviour, Q ^ T. In a model 
of Id conductors with dilute weak-links, a linear metal- 
lic thermopower would be expected. The Lorenz number 
provides a further difference between these two models. 
At the Anderson transition the Lorenz number is pre- 
dicted to be supressed ||l9] below the Fermi-liquid value 



Lq by roughly 2/3, whereas our results show an enhanced 
Lorenz number for Id thermal transport through dilute 
impurities. It would be most interesting to measure ther- 
mal conductivity in conducting polymer samples which 
exhibit power law electrical conductivities. 
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